Design of effective kernels for spectroscopy and molecular transport: time-dependent 

current-density- functional theory 
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Time-dependent current-density-functional theory (TDCDFT) provides an in principle exact 
scheme to calculate efficiently response functions for a very broad range of applications. How- 
ever, the lack of approximations valid for a range of parameters met in experimental conditions 
has so far delayed its extensive use in inhomogeneous systems. On the other side, in many-body 
perturbation theory (MBPT) accurate approximations are available, but at a price of a higher com- 
putational cost. In the present work the possibility of combining the advantages of both approaches 
is exploited. In this way an exact equation for the exchange-correlation kernel of TDCDFT is 
obtained, which opens the way for a systematic improvement of the approximations adopted in 
practical applications. Finally, an approximate kernel for an efficient calculation of spectra of solids 
and molecular conductances is suggested and its validity discussed. 



PACS numbers: 71.10.- 



71.45.Gm, 78.20.-e 



The theoretical description of the response of an elec- 
tronic system to a time-dependent perturbation is a key 
problem for many areas of physics and chemistry. Most 
spectroscopic experiments probe the elementary excita- 
tions of an electronic system through its linear response 
to an external electromagnetic field. Their theoretical in- 
terpretation is of primary interest for technological appli- 
cations in condensed-matter physics, nanosciences, pho- 
tochemistry or biophysics. Similarly, response functions 
are essential, among many possible applications, also for 
the determination of the electrical conductivity or other 
transport coefficients in molecular electronic devices. 

Therefore, one would like to devise a reduced the- 
oretical framework that is, at the same time, reliable 
and efficient, by calculating only the information needed 
to interpret and predict specific experimental measure- 
ments. Two prominent examples of such reduced ap- 
proaches are many-body perturbation theory (MBPT) 
and density-functional-based theories. Key variables of 
the former are one- and two-particle Green's functions, 
G(l,2) and G2(l,2, 3,4) (1 is a shorthand notation for 
space, time and spin indices V\,ti,a\). Methods based 
on the Green's-function formalism reduce the complex- 
ity of the many-body wavefunction into the propaga- 
tion and the interaction of renormalized quasiparticles. 
Their intuitive, direct contact with the initial problem 
of real interacting electrons makes it rather easy to in- 
troduce working approximations. A remarkably suc- 
cessful example is the solution of the Bethe-Salpeter 
equation (BSE) for the two-particle correlation function 
£(1,2,3,4) = -G 2 (l,2,3,4) + G(1,3)G(2,4) which has 
led to an important breakthrough by permitting, for in- 
stance, an accurate calculation of electronic spectra of 
solids and nanosystems^. On the other hand, practi- 
cal calculations at this level remain very demanding, 
even for nowadays' computers. An alternative pathway 



is instead based on the extension of density-functional 
theory (DFT) to scalar time-dependent external poten- 
tials, V ex t(r, t), as in time-dependent density- functional 
theory (TDDFT) 2 , or to time-dependent vector poten- 
tials, A ext (r, t), as in time-dependent current-density- 
functional theory (TDCDFT) 3 ' 4 . When one needs only 
charge or current densities, p(r,t) and j(r, t), these 
density-based methods identify the minimum content of 
information that one has to calculate in order to provide 
the searched answers. In the Kohn-Sham (KS) scheme^, 
the many-body problem is reformulated very efficiently 
into a set of self-consistent non-interacting one-particle 
equations. For this reason, in the KS scheme the solu- 
tion of the full many-body problem is made simple and its 
computational cost is very convenient. The main draw- 
back is that it is generally very difficult to improve upon 
the simplest local-density approximations (LDA)^ for the 
exact density functionals of the formal theory. In fact, 
after the first promising results^, the development of TD- 
CDFT in this field has been delayed by the lack of ade- 
quate approximations to the tensor exchange-correlation 
(xc) kernel f xc , beyond the local functional in the cur- 
rent density of Vignale and Kohn (VK)^£, which is un- 
fortunately valid only in a range of parameters that is 
often not met in experiments performed on inhomoge- 
neous systems^. 

The present work aims at overcoming the limitations of 
the VK functional, by applying to TDCDFT the emerg- 
ing successful strategy of combination of MBPT and 
density-functional approaches, in order to profit from 
the complementary advantages of bothlSr— . In particu- 
lar, we will derive an exact relation linking the unknown 
exchange-correlation tensor kernel f xc to quantities that 
can be in principle accurately calculated in MBPT. More- 
over, we will show how it is possible to introduce suitable 
approximations to this exact relation opening the way to 
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a broad class of applications and to a systematic way to 
improve the approximations adopted in TDCDFT calcu- 
lations. 

The Bethe-Salpeter equation for the irreducible po- 
larization function 1/(1,2,3,4) reads 15 (throughout the 
paper, integrals and sums are always done on repeated 
indices and atomic units are adopted): 

1(1,2,3,4) = £,,(1,2,3,4)+ 

L (1,2,5,6)S(5,6,7,8)L(7,8,3,4). (1) 

Here L (l,2, 3, 4) = -iG(l,3)G(4,2) is the two- 
particle correlation function for independent particles 
and 3(5,6,7,8) = i5E(5, 6)/<5G(7, 8) is the BSE ker- 
nel, which e.g. accounts for excitonic effects in op- 
tical spectra. In standard BSE implementations the 
GW approximation-^ for the self-energy E is adopted, 
where E is evaluated as a product of the one-electron 
Green's function G and the screened Coulomb interac- 
tion W. S is usually approximated as 2(5,6,7,8) = 
— W(5, 6)<5(5, 7)6(6, 8). This amounts to neglecting the 
term iGSW/SG, which contains information about the 
change of the screening in the excitation and is consid- 
ered to be small. Moreover, for W one generally con- 
siders only a static screening of the Coulomb interac- 
tion v and Lq is built with GW quasiparticle (QP) en- 
ergies and KS wavefunctions. L is then linked to the 
correlation function L by a Dyson equation: L = L + 
LvL. Whereas the quantities of spectroscopic interest are 
for instance the two-point density-density and current- 
current response functions, x pp (l,2) = S p(l) / SV ex t{2) 
and x(l,2) = Sj(l)/6A ext (2), the BSE is an intrinsi- 
cally four-point equation. In fact, in the BSE scheme, 
these two-point response functions can be obtained only 
as contractions of the four-point correlation function L, 
which has to be calculated in a first step. In many situ- 
ations, as the ones we are interested in here, this clearly 
reveals to be a computational waste that one would like 
to avoid. 

In TDCDFT the linear response of the current j to an 
external vector potential A ext is^: 



$?«(!) = ^,^(1,2)^(2), 



(2) 



where the Kohn-Sham vector potential A s is the sum 
of the external, Hartree and exchange-correlation poten- 
tials: A s (l) = Aexi(l) + A ff (l) + A xc (l), and Xs is 
the Kohn-Sham current-current response function. Simi- 
larly, the linear variation of the current can be calculated 
through the knowledge of the irreducible current-current 
response function \: 

SjM) = \xap{h 2)[8A ext ^(2) + 6A Ht0 (2)]. (3) 

Combining these two definitions, one immediately gets 
to a Dyson equation linking Xs with \' 

Xa/8(l,2) = X«,a/?(l,2)+ 

Xs,aA(l,3)/, c , AK (3,4)^(4,2), (4) 



where the exchange-correlation tensor kernel 

8A XCj0l (l) 



fxc,afl(l, 2) 



tys(2) 



(5) 



has been introduced. Once thanks to (j4} the irreducible 
X has been calculated, the (reducible) response function 
X (hence the spectra) can be obtained through: 



Xa/3(ri,r 2 ,w) = x a/3 (ri,r 2 ,cj)+ 



Xax(ri,r 3 ,o;)- 



l r 3 - r 4 



rV4 K X K /j(r4,r 2 ,w). (6) 



The quality of the approximation adopted for the xc ker- 
nel ([5]) is hence fundamental for the accuracy of the final 
results. 

As in TDDFT, also in TDCDFT the effect of the 
xc kernel on the spectra calculated from the indepen- 
dent KS-particle response is twofold. The Kohn-Sham 
eigenvalues are known to underestimate the quasiparti- 
cle band gap of insulating systems due to the derivative 
discontinuity of the DFT xc potential^. Hence the xc 
kernel has first to provide a consistent band gap open- 
ing. And, second, as the BSE kernel S, it has to de- 
scribe electron-hole interactions. So, following Refsi 12 ' 13 , 
in order to make explicit these two aspects, here we set 

f X c,aj3 



J xc,aj3 



Sap. where fxLf) has the task of 



(2) 



for 



overcoming the KS band gap problem, while /, 
instance, accounts for excitonic effects in optical spec- 
tra, or dynamical corrections to the Landauer formula 
for the electronic conductance in quantum transport 18 . 
Formally, we will now introduce a contraction operator 



U af} Lo(l,l',2,2' 

1 1 kv, 



2i 2i 



Vv)(V 2 , -V 2 jLo(l,l',2,2') 



= l+,2'=2+ 



(7) 



in such a way that xo 

X0M 1 , 2) = p(l)5(l, 2)S aP + n Q/3 L (l, I', 2, 2'). (8) 

Xo is built with QP ingredients instead of KS ones. In 
this way one has: 



[XJ 1 ]^ (1,2) - [xo 1 ]^ (1,2) = /£L,(1,2), 
and, then: 

(1=2) = £^(1,2). 



r 1 



a/3 



(9) 



(10) 



In many semiconductors, the difference between KS and 
QP can be accounted for by using a scissor operator that 
shifts rigidly upwards the eigenvalues of the conduction 
with respect to valence bands. The use of a scissor opera- 
tor is also a common practice in TDCDFT-^. The main 
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point of interest in our discussion here is hence about 

f (2) 

Eq. (TO) leads to: 



^(l,2) = xo, Q /3(l,2)+ 

Xo, Q A(l,3)/i 2 c | AK (3,4fx K/3 (4 ! 2). (11) 

By definition, both TDCDFT and BSE yield the exact 
two-point response function Xa/3(1, 2) = p(l)S(l, 2)5 a /3 + 
n a pL(l, 1', 2, 2'). Therefore, applying the contraction 
n Q( g to Eq. (fTJ) and comparing the result with Eq. (fTT]) . 
one obtains: 



*o,oa(1, 3)^(3, 4)^(4, 2) = 

□ Q/3 {l (1,1',3,4)S(3,4,5,6)L(5,6,2,2')}, ( 

m 



12) 



which can be solved for the TDCDFT kernel f^. c a a giv- 
ing: 



/ic!^(l. 2) = [Xo '] aA (l, 3) n AK [£ (3, 3', 4, 5) 



3(4,5,6,7)1(6,7,8,8') [x 



(13) 



Eq. ([12]) is a generalized Sham-Schliiter equatio n 17 ' 20 , 
which relates in an exact manner TDCDFT quantities 
with MBPT ones, opening the way to possible system- 
atic improvements in the design of new approximations 
for the xc kernel of TDCDFT. Its main worth is that it 
allows avoiding to approximate directly the kernel f xc , 
which can be a difficult problem. The design of opera- 
tive approximations is instead simpler in the context of 
MBPT. Then, thanks to Eq. (fT3"|) . working approxima- 
tions of MBPT can be mapped into the more efficient 
TDCDFT scheme, where one would prefer to solve the 
equations. In particular, when approximating Eq. (p~3|) . 
no assumptions of locality in the current density of the xc 
functional are explicitly needed, leading to approximate 
kernels that can be employed also in the range of pa- 
rameters where the VK functional is formally not valid, 
namely for ground-state densities and induced current 
densities not slowly varying in space, and in the region 
below the particle-hole continuum of the homogeneous 
electron gas. 

Within the TDDFT framework, such a mapping strat- 
egy has already demonstrated to be a successful approach 
and has led to the introduction of an xc kernel, known as 
Nanoquanta kerneU°- ~ 14 ' 20 , which has shown to provide 
the same level of accuracy as the BSE in a wide range of 
spectroscopy applications, from solids to finite molecular 
chains 21 ' 22 . Therefore, by discussing a first practical ap- 
plication of the exact Eq. (|13|) . here we will show that 
a similar approach for the design of new approximations 
to the tensor xc kernel of TDCDFT is also possible. In 
fact, one expects to find a similar level of accuracy also 
for current-current response functions. By a first-order 
linearization of Eq. (I13|) , where for the various response 



functions one uses the independent-particle versions, and 
taking for 3 the statically screened W, as usually done 
in BSE, one obtains: 



More explicitly: 



/Sa /3 (l ! 2) = -[Xo 1 ] aA (l ; 3) 



[Xo 1 ] QA (1, 3)n AK [L (3, 3', 4, 5)W(4, 5) 
Lo(4,5,6,6')][x^]^(6,2). (14) 



^(V 3a -V 3 <)£ (3,3',4,5)^(4,5) 



lim 

3'->-3 

6'-+6 



"(Ve, -V 6; )Lo(4,5,6,6') 
2,1 



[Xo'UCM, (15) 



where the presence of three-point current-density and 
density-current response functions, x| p (3;4, 5) and 
Xoj(4, 5; 6), becomes apparent. Hence the kernel (|15[) can 



be also rewritten in a compact way as: 



(2) _ 



-Xo^ljpWx^pjXo 1 - 



(16) 



This is a new approximation of the TDCDFT f xc kernel 
that has to be understood as an orbital functional, hence 
an implicit functional of the current density. In this sense 
it is more flexible than the VK functional, which instead 
is an explicit functional of the current density. 

The spatial derivatives that appear in Eq. ([15)) don't 
modify the structure of the poles of the response func- 
tions in the frequency domain. Therefore, this approx- 
imation of the TDCDFT tensor xc kernel can benefit 
from the same cancellation of poles and zeroes of the re- 
sponse functions that has been shown to be essential for 
the Nanoquanta kernel of TDDFT 2 ^ (provided that the 
response functions entering Eq. ([TBI are built with QP 
energies). Since any scalar potential with a gauge trans- 
formation can be represented by a longitudinal vector 
potential, Eq. (TTBl can be thought as a generalization 
of the Nanoquanta TDDFT f xc kernel to the calculation 
of the response to any kind of time-dependent external 
vector potential. 

"(2) 

In general the relation between the tensor f xc ' TD- 
CDFT kernel and the scalar TDDFT kernel is rather 
involved 2 ^. Its first-order linearization in fxc read 



,25. 



(2) 

XC 



Xo.LVxo/i^XoVXoi 



pp 



(17) 



This approximation is consistent with the linearization 
of Eq. ([13")) that leads to (fTBf . Hence, by inserting the 
TDCDFT kernel ^ in Eq. (JT7|) and using the 



density continuity equation 2 ^, the TDDFT f xc ' 
becomes: 



kernel 



which is, in a compact form, the Nanoquanta TDDFT 
kernel ^ 14 ' 20 . This equivalence, to the first order in 
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fxc , between the tensor kernel (fTB"|) and the scalar ker- 
nel (|18|) supports the validity of the new approximation 
for the TDCDFT kernel ^T5J), thanks to the excellent re- 
sults found using the Nanoquanta TDDFT kernel ([TBF 1 . 
Moreover, this represents a further alternative derivation 
of the expression ([T8l) of the Nanoquanta kernel that ex- 
ploits the possibility of mapping approximations devel- 
oped in the framework of TDCDFT into scalar TDDFT 
kernelo 24 - 25 ! 27 - 28 . 

In situations where the screening of the Coulomb in- 
teraction is ineffective, W can be approximated with v 
and the GW approximation reduces to Hartree-Fock. In 
this case, (|15l) reduces to an exact-exchange approxima- 
tion for f xc (which has been already worked out for the 
homogeneous electron gas in RefJ^). On the other side, 
it has been demonstrated within the TDDFT framework 
that the lack of screening of the long-range contribution 
of the Coulomb interaction in the kernel implicitly over- 
estimates both QP band gaps and excitonic effects, lead- 
ing to pathologies in optical spectra of semiconductors 2 ^. 
These pathologies can be cured by taking into account 
the screening of the Coulomb interaction, as done in Eq. 

G5). 

For optical spectra of solids, relevant is the long- 
wavelength limit q — > of the G = G' = element of the 

"(2) 

xc kernel written in reciprocal space: f xc a/3 (q + G, q + 
G',ui), where G and G' are reciprocal lattice vectors. 
In this limit the f x 2 ) kernel (|T5j) in insulators becomes^ 
a/uj 2 . Moreover, the static W is proportional to l/e^, 
where is the static dielectric constant. This implies 
that for the q — > limit also a is proportional to l/foo, 
suggesting that in isotropic systems a could be used as a 
fitting parameter for the calculation of optical spectra in 
solids 10,31 . In this limit, the tensor kernel becomes com- 
pletely local, contrary to the scalar TDDFT kernel that 
is ultranonlocal 32 and has a 1/q 2 asymptotic behavior in 
the long-wavelength limi t 10 i 33 . 

Here it is interesting also to note that, even though for 



the BSE kernel a static approximation for W is adopted, 
the resulting f xc kernel in (1151) in general is naturally 
frequency dependent^. This is a consequence of the 
conversion of the spatial nonlocality into a frequency 
dependence 2 ^, associated to the reduction of the number 
of degrees of freedom, when one passes from the four- 
point L to the two-point x an d from the four-point 5 to 
the two-point f xc . 

Since the validity of the new f xc kernel (fT5|) is not 
confined to the weakly inhomogeneous limit, it could be 
for instance used also to study dynamical corrections in 
the weak bias limit of the Landauer formula in molecular 
transport for nanoscale junctions^. In fact, so far the 
calculated corrections have been based on the VK func- 
tional, which is strictly valid only for slowly- varying den- 
sities in a high-frequency regime. Therefore an accurate 
estimate of these effects is still under debate 3 ^. On the 
other side, the renormalization of molecular electronic 
KS levels at metal-molecule interfaces^ in the present 
context is accounted for by the term fxc . 

In conclusion, in the present work we have derived an 
exact equation for the xc kernel f xc of TDCDFT, which 
allows one to map well established working approxima- 
tion of MBPT into the more efficient TDCDFT scheme, 
Due the lack of an appropriate approximate xc kernel, 
so far TDCDFT could not be used for calculating op- 
tical spectra in good agreement with experiments^^. 
So, in analogy with the successful results obtained in the 
TDDFT case 2 ^, we have suggested a practical approxi- 
mation for f xc and discussed its validity for the calcula- 
tion of electronic spectra of solids and nanosystems and 
dc conductances in molecular devices. 
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